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Abstract 

By dividing the redshift region under consideration into two bins, we estimate the limitations 
in determining the equation of state Wd e of dark energy at high redshift from current and future 
observational data including supernovae, baryon acoustic oscillation and observational Hubble data. 
It is found that the constraints of Wd e from current data are weak ( 2a(wd e ) ~ 1 ) even beyond 
redshift z = 0.3. For simulated future observational data, it is shown that ~ 2300 supernovae 
data from a SANPdike JDEM survey give 2a(wd e ) ~ 1 beyond z = 0.6. We consider the effects 
of the divided point, the number of supernovae data and the error in the distance modulus on 
constraining Wde at high redshift. It is shown that the increase of number of supernovae data 
seems not efficient on improving the constraints of u>d e at high redshift, while the improvement of 
observational error in distance modulus seems to be a better way. 
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I. INTRODUCTION 



The current expansion of the universe is found to be in accelerating jl , 2 ] and it is believed 
that the dark energy (DE), whose equation of state Wd e = Pdel Pde is less than —1/3, plays 
the role to drive the accelerated expansion of the universe. Plenties of DE models have been 
proposed 3|-[9| but the origin of DE is still unknown. One has only known a little about 
DE from observational data. Several parametrization methods of Wd e have been proposed 
to fit with observations to get constraints of Wd e - Usually in parametrization % is assumed 



to be a simple function of redshift z, e.g., the CPL parametrization [lO, 11|: Wd e (z) 
Wo + w a z/(l + z), and the parametrization of redshift expansion [l^Q: Wd e (z) = w + w z z. 
It has been found that Wd e is very close to —1 at low redshift and is varying very slowly (if it 
is dynamical). Of course the fitting results depend on the parametrization forms of Wd e {z). 
And an assumed form of Wd e {z) might be only suitable to describe one type of w de - For 
example, the CPL parametrization cannot describe the equation of state with oscillations 
or rapid transitions. In general, one parametrization form of Wd e (z) could only approximate 
the real Wd e well in a very limited region of redshift. 



Several model- in dep endent methods have also been proposed to reconstruct 



observational data 
spline method 



14 



17], e.g., uncorre 



201 ]. linear-spline method 



ated band-power estimates (UBE) [18|, [l9| , cubic 



z) from 



2 1| .wavelet approach 



22 1 , and so on. Most of these 



methods relate to a piecewise Wd e (z): 



w de (z) = Wi(z) , Zi < z < Z i+ i 



(1) 



where zq = and Wi(z) is a simple function of redshift z in the i bin. It is not surprisin 



that in higher redshift bins, constraints of w^e from observational data will be weaker [23]. 
There are mainly two reasons: 1) In higher redshift bins the number of currently available 
data points will be much less. At present, the most data points to constrain DE are from 
Type la supernovae (Snla). There are much less Snla data points in higher redshift bins, as 
shown in Table I. Moreover, w de in one bin can also be constrained by data points of other 
bins that have higher redshift, i.e., w de in the lowest redshift bin can be constrained by all 
data points. 2) In higher redshift bins, the fractional energy density of DE Qde = Pde/3H 2 
is always be smaller than the one of dark matter, so the role played by DE in the luminosity 
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TABLE I: The redshift distribution of 472 supernovae data from 3 years Supernova Legacy Survey. 
The redshifts in the first row are the upper limits of each bin. 



distance 



Di(z) = {l + z) / dx/H(x) 
Jo 



(2) 



will be weaker. This can be seen in a simple calculation as follows. In a flat FRW universe 
the Fridemann equation reads: 



where 



# 2 0) = g[Pmo(l + xf + p de (x)}, 



Pde{x) = p de0 e Jo *+» 



(3) 



(4) 



Supposing there is a tiny variation 5wi of Wd e in the i th bin and no variations of Wd e in all 
other bins, from Eq.([2|) and ([3]) one has 



0, 



< x < z,; 



SH 2 (x) = < 



p de In {j^r)5Wi , Zi < X < Z i+ i 



■ 1 + 2 

Pcfe ln(i±gi)^ 



(5) 



Zi+i < x 



and correspondingly 



5D l (z) 



ST In (ig)dz + !' ^ de In , 



l + Zi. 



l+2i 



■)da; 



< z < z.i 

Z{ < Z < Zi+i 

(6) 



We can see that 5wi always appears together with fide in 8D[(z). Since VL de is much 
smaller at higher redshift, SDi(z) is also a smaller change in higher redshift bins with the 
same 5wi, i.e., 5Di(z) will be less sensitive to 5wi. It indicates that when the redshift in a bin 
is larger (i is larger), to lead to the same variation of Di the change of Wi must be larger. In 
other words, the constraint on Wd e from the luminosity distance is weaker in higher redshift 
bins than in lower redshift ones. 
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As shown in Eq.([6]), 5Di in the i th bin is an integration, which is mainly effected by 
evolving of fide with z and the width of this bin. At high redshift bins Vide will be small and 

n de ( Zi < z < z l+1 ) ~ (1 + zf Wi . (7) 

When Wi is larger and the bin is wider, 5Di in the i th bin will be larger. There is a situation 
which needs to be mentioned that at a high redshift bin if is much smaller than —1, the 
decrease of fide will be so rapid that 5Di keeps extremely small even the bin is very wide 
and Swi is very large. In this case, the likelihood of Wi will be very flat in a large region 
of small Wi as shown in Fig. [2] in the next section, since Wi in this large region almost has 
the same D[. Thus the downward errors of Wi at high redshift bins are always extremely 
large. In general, the constraint of Wd e in the last redshift bin is always the weakest and the 
downward error of Wd e is always very big. 

In this paper, we would like to estimate the limitations in determining Wd e at high redshift 
from present and future observational data. Constraints of Wd e will be obtained by using the 
UBE method with 2 redshift bins and errors of Wd e in the second bin indicate the limitations 
in determining Wd e at high redshift. The effects of divided positions of two redshift bins, 
the number of Snla data and observational error of Snla data in constraining Wd e at high 
redshift will be studied. 



II. CONSTRAINTS OF w de FROM PRESENT DATA 



We will use the UB 
3 years data (SNLS3) 



method to get constraints of Wd e from 472 Supernova Legacy Survey 



24j, Baryon Acoustic Oscillation (BAO) Distance Parameter A 25 ] 



2J28 



. In this case, redshift for all 



and 12 Observational Hubble Data with z < 1.4 from 
those data points are smaller than 1.4 and one needs not to involve Wd e in higher redshift 
region with z > 1.4. In Table [Til 12 Hubble data are listed. And the BAO Distance 
Parameter A is defined as: 



, = fi V 2£ -V3 (0 . 35) (_i_^^.) 



2/3 



(8) 



where E(z) = H(z)/H(0) and Q m0 = p rn0 /3H 2 (0). The redshift region of [0,1.4] will be 
divided into two bins and 

W\ , < z < Z\ 
w 2 , Z\ < z < 1.4 



Wde(Z) 



(9) 
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TABLE II: 12 observational Hubble data with their redshifts z < 1.4 from 
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where w\ and W2 are two constants. The divided position z\ will be set as 0.3, 0.35, 0.4 
and 0.45, respectively, as shown in Fig. [TJ The best-fitted parameters and their 68.3%, 
95.4% confidence level (C.L.) errors are obtained by using the Markov chain Monte Carlo 
(MCMC) method, which are shown in Table IIII1 The errors of u>2 can be used to estimate 
the limitations in determining Wde beyond z — Z\. In the numerical calculation we have set 
the prior with w 2 > —20, or u> 2 would go to a very large minus value in MCMC procedure 
and the downward error of u>2 would be extremely large, as discussed in the introduction. 
The figure of merit (FoM) defined as {29], I30I ]: 



FoM = [detC(u>i,u> 2 )] 



-1/2 



(10) 



has also been calculated for each model, where C(wq, w\) is the covariance matrix of w\ and 
W2 after marginalizing out all other parameters of the model. Moreover, it is found that the 
correlation between w\ and W2 is very small, i.e., the errors of w\ and W2 obtained can be 
treated as independent with each other. 
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TABLE III: The best-fitted values and their 68.3% and 95.4% C.L. errors of the UBE models from 
present observational data. In the UBE models, the redshift region (0, 1.4) is divided into two bins 
and the divided positions are z\ = 0.3, 0.35, 0.4, 0.45, respectively. The figure of merit in each case 
is also calculated. 



From the above we may conclude as follows. 

1. When the redshift region of the second bin becomes smaller, the constraints of wj 2 
then become much weaker. As shown in Fig. [fl the errors of W2 increase very rapidly with 
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FIG. 1: The best-fitted w de {z) and its 68.3% and 95.4% C.L. errors of the UBE models from the 
present observational data. 

the divided point Z\. It is shown in Table [TTT1 that the errors of Wd e beyond Z\ = 0.3 is large 
(the width of uncertainty of W2 at 95.4% C.L. is 1.93, i.e., 2a (w 2) ~ 1), compared with the 
ones with z < 0.3. Furthermore, if one sets up Z\ > 0.45, the constraints of Wd e beyond z 1 
will be extremely large (its lower limit at 95.4% C.L. is smaller than -20). Therefore it is 
very hard to determine w^e at high redshift from current observational data. 

2. At high redshift, the likelihood of Wd e is flat on its downward side. As shown in Fig. [2j 
for larger z\ the likelihood of W2 is flatter on the downward side, and thus the downward 
errors of W2 are much larger, as shown in Fig. [TJ 

3. The effect of the divided point Z\ on Wi, i.e., Wd e in low redshift region, is very small. 
With very different z\ there is a little difference between the errors of corresponding w\, as 
shown in Table III II and Fig. [TJ It indicates that the number of present data at low redshift 
is sufficient and constraints on Wd e in this region is strong. 

4. FoM can be used to estimate the goodness of constraints of Wd e at high redshift from 
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FIG. 2: The likelihood of w 2 of the UBE models with the divided points z x = 0.3, 0.35, 0.4, 0.45, 
respectively. 
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TABLE IV: The redshift distribution of 1998 SN data for a SNAP-like JDEM survey and 300 SN 



data from the NSNF 



31 



321 ] . The redshifts in the first row are the upper limits of each bin. 



observations. FoM from Eq.f JTOj) is proportional to the inverse area of the error ellipse in 
the Wi ~ u>2 plane. Since errors of w\ change very slowly with z\, differences between FoM 
with different z\ are mainly from differences between the errors of w%. FoM can be used to 
estimate the effect of Z\ on constraints of w 2 , i.e., the goodness of constraints of Wd e beyond 
z = z\ . As shown in Table IHIt FoM increases rapidly with the decrease of z\ . 



III. CONSTRAINTS OF w de FROM FUTURE DATA 

To estimate the constraints of Wd e at high redshift from future observations, we will fit 
the 2-binned UBE model shown in Eq.([9]) with 2298 simulated supernovae data (denoted as 

7 



elementary data set in this paper), which contain 1998 supernovae data with 0.1 < z < 1.7 
rom a SNAP-like JDEM survey and 300 supernovae data with z < 0.1 from the NSNF 



31 



32| . The redshift distribution of the mock supernovae data is shown in Table IIVI To 



331 ]. the data of the BAO Distance Parameter 



alleviate the degeneracy between fl m0 and w& 
A 25] will also be included. To simulate the mock supernovae data, we assume the fiducial 
model as w ^(z) = — 1. The error in the distance modulus for each supernovae will be set 



as 5 = 0.13 [17|, except in the subsection C. In this section we will study the effect of the 
position of z\, the number of supernovae data and the value of 6 in constraining Wd e at high 
redshift. 



A. The effects of the divided position of bins 
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TABLE V: The best-fitted values and their 68.3% and 95.4% C.L. errors of the UBE models from 
the 2298 simulated supernovae data. In the UBE models, the redshift region (0, 1.7) is divided into 
two bins and the divided positions are z\ = 0.4,0.6,0.8,1.0, respectively. The figure of merit in 
each case is also included. 

To estimate the effect of the divided position z\ on errors of Wd e , we have set z\ as 0.4, 
0.6, 0.8 and 1.0, respectively, as shown in Table IVl and Fig. [3j As expected, errors of u^ e at 
high redshift (i.e., ^2) increase rapidly with z±. One can see that constraints of Wd e is weak 
beyond z = 0.6: in the case of z\ = 0.6, the width of error region of u>2 at 95.4% C.L. is 
1.76 (i.e., 2a(w2) ~ 1), which is similar to the constraints of Wde beyond z = 0.3 from the 
current data (in the previous section). It also shows that constraints of Wd e beyond z = 1 
becomes extremely weak. On the other hand, the constraint on Wd e at low redshift (u>i) will 
be much better than that from the present data. Finally one can see that the value of Xmin 
changes very slowly with z\. 
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FIG. 3: Plots of the best-fitted Wd e (z) and their 68.3% and 95.4% C.L. errors shown in Table El 
B. The effects of the number of supernovae data 
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TABLE VI: The best-fitted values and their 68.3% and 95.4% C.L. errors of the UBE models from 
simulated supernovae data. The numbers of supernovae data are respectively 2, 3, 4, and 5 times of 
the 2298 elementary data, respectively. In the UBE models, the redshift region (0, 1.7) is divided 
into two bins and the divided position z\ is set to 1.0. The figure of merit for each case is also 
calculated. 

In the subsection A it is shown that ~ 2300 mock future supernovae data still give poor 
constraints of Wd e beyond z = 0.6. Here we would like to see the effects if we increase the 
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FIG. 4: Plots of the best-fitted Wde(z) and their 68.3% and 95.4% C.L. errors shown in Table IVH 

number of the mock supernovae data, with 2, 3, 4 and 5 times of the 2298 elementary data 
introduced in the subsection A. The proportional distribution of redshift in these multiple 
datasets is the same as the elementary data set, as shown in Table HVl The divided position 
of two bins z\ is set to 1.0 . The results are shown in Table I VI I and Fig. HI It shows that 
the improvement of FoM (i.e., constraints of Wde beyond z — 1) with increase of the number 
of data points is not very efficient. Constraints of w^e beyond z = 1 from the data with 5 
times number of the elementary supernovae data set are still weak (2a (W2) ~ 1 with 5 times 
number of data). 



C. The effects of the error in the distance modulus 



One can see that increasing the number of supernovae data seems not efficient enough to 
improve constraints of Wd e at high redshift. However, it is expected that with the increase 
of number of data, the statistical errors will be improved and the systematic errors should 



also be improved more and more by future observations 34J. To estimate effects of the error 



S, here we set 8 = 0.13, 0.1, 0.05, 0.02, respectively. The number of the supernovae used here 
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TABLE VII: The best-fitted values and their 68.3% and 95.4% C.L. errors of the UBE models with 
errors in distance modulus 5 = 0.13,0.10,0.05,0.02, respectively. The number of supernovae used 
is 3 times of 2298 elementary supernovae data. In the UBE model, the redshift region (0, 1.7) is 
divided into two bins and the divided positions z\ is set to 1.0. The figure of merit for each case 
is also included. 
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FIG. 5: Plots of the best-fitted w de (z) and their 68.3% and 95.4% C.L. errors shown in Table EH 

will be 3 times of the elementary data and z\ = 1.0 is set. The results show in Table EIH that 
FoM increases very rapidly with the decrease of S. For example, as S = 0.02, the constraints 
of Wde in z G (1, 1.7) is even better than the present constraints of w^e at low redshift bin 
obtained in previous section. Moreover, Wd e at the low redshift region can be definitively 
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determined with 5 = 0.02, as shown in Table IVlTI 



IV. SUMMARY 

We have used the UBE method to estimate the limitations in determining Wde at high 
redshift from current and future mock supernovae data. The redshift region of the data 
points is divided into two bins with the divided position z\. Best-fitted wi, W2 and their 
errors at 68.3%, 95.4% CL. are obtained by using the MCMC method. FoM for each case 
is also calculated, which can be used to estimate goodness of constraints on 102. Constraints 
on W2 represent the ability of observational data in determining Wde beyond z — Z\. 

For the present observational data set, 472 Supernovae Legacy Survey 3 years data, BAO 
Distance Parameter A, and 12 Observational Hubble data with z < 1.4 are adopted. It 
has been shown that the errors of w<i increase rapidly with the divided position Z\. For the 
current observational data, the constraints of Wd e beyond roughly z = 0.3 are still very weak. 
The downward errors of Wd e at high redshift are extremely large, since the likelihood of Wd e 
at high redshift is very flat on the downward side. 

For the future data, we have simulated supernovae data from a SNAP-like JDEM survey. 
The effects of divided position Z\, the number of supernovae data points and the error in the 
distance modulus are analyzed. It was shown that the constraints of Wd e from the mock 2298 
supernovae data are extremely weak beyond roughly z — 1. The increase of the number of 
supernovae data seems not as efficient as the decrease of the error of distance modulus in 
improving constraints of Wd e at high redshift. 

Nowadays, many projects for the observations of DE are in progress or in planning, 
it is expected that a great deal of data points would be released. The observation from 
supernovae is still the most important one to reveal the behavior of DE. To make sure 
whether Wd e is dynamical or not, a narrow constraint of Wd e at high redshift is required. We 
need much more data points, and the improvement of accuracy of the observations seems 
more necessary. It is shown in Table IVHI that ~ 6900 Snla data with their observational 
errors o = 0.02 will constrain uncertainties of Wd e beyond z = 1 within 0.1 at 68.3% CL. 
and within 0.2 at 95.4% CL. . Note that all observations used and analyzed in this paper 
relate to the distance measurement. Obviously it is believed that more observations of other 
types for DE combining with distance measurements should improve the constraints of Wd e 
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at high redshift. 

Acknowledgements: 

This work was supported in part by the National Natural Science Foundation of China 
(No.10821504, No.10975168, No.11035008 and No.11147186), in part by the Ministry of Sci- 
ence and Technology of China under Grant No. 2010CB833004, and a grant from Hanzhou 
Normal University. 



A. G. Riess et al. [Supernova Search Team Collaboration], Astron. J. 116 (1998) 1009 
astro- ph/9805201] . 

S. Perlmutter et al. [Supernova Cosmology Project Collaboration], Astrophys. J. 517 (1999) 
565 [astro-ph/9812133] . 



P. J. E. Peebles and B. Ratra, Rev. Mod. Phys. 75 (2003) 559 [astro-ph/0207347] . 

E. J. Copeland, M. Sami and S. Tsujikawa, Int. J. Mod. Phys. D 15 (2006) 1753 

|hep-th/0603057] . 

P. J. Steinhardt, L. -M. Wang and I. Zlatev, Phys. Rev. D 59 (1999) 123504 
[astro-ph/9812313] , 

S. Capozziello, S. Carloni and A. Troisi, Recent Res. Dev. Astron. Astrophys. 1 (2003) 625 
astro- ph/0303041] . 

R. -G. Cai, Phys. Lett. B 657 (2007) 228 |arXiv:0707.4049l [hep-th]]. 

R. -G. Cai and Q. Su, Phys. Rev. D 81 (2010) 103514 |arXiv:0912. 19431 [astro-ph.CO]]. 

R. -G. Cai, Z. -L. Tuo, H. -B. Zhang and Q. Su, Phys. Rev. D 84 (2011) 123501 

[arXiv: 101 1.32121 [astro-ph.CO]]. 



M. Chevallier and D. Polarski, Int. J. Mod. Phys. D 10 (2001) 213 (grjqc/0009008]. 
E. V. Linder, Phys. Rev. Lett. 90 (2003) 091301 |astro-ph/0208512] . 
A. R. Cooray and D. Huterer, Astrophys. J. 513 (1999) L95 [astro-ph/9901097] , 
D. Huterer and M. S. Turner, Phys. Rev. D 64 (2001) 123527 [astro-ph/0012510] . 
D. Huterer and G. Starkman, Phys. Rev. Lett. 90 (2003) 031301 |astro-ph/0207517] . 



D. Huterer and A. Cooray, Phys. Rev. D 71 (2005) 023506 |astro-ph/0404062| . 

Y. Wang, Phys. Rev. D 80 (2009) 123525 |arXiv:0910.2492l [astro-ph.CO]]. 

T. Holsclaw, U. Alam, B. Sanso, H. Lee, K. Heitmann, S. Habib and D. Higdon, Phys. Rev. 



13 



D 82 (2010) 103502 [arXiv: 1009.54431 [astro-ph.CO]]. 
i] S. Sullivan, A. Cooray and D. E. Holz, JCAP 0709 (2007) 004 |arXiv:0706.3730l [astro-ph]]. 
3] G. -B. Zhao, D. Huterer and X. Zhang, Phys. Rev. D 77 (2008) 121302 |arXiv:0712.2277l 
[astro-ph]]. 

3] P. Serra, A. Cooray, D. E. Holz, A. Melchiorri, S. Pandolfi and D. Sarkar, Phys. Rev. D 80 

(2009) 121302 |arXiv:0908.3T86l [astro-ph.CO]]. 
1] R. -G. Cai, Q. Su and H. -B. Zhang, JCAP 1004 (2010) 012 [arXiv:1001.2207l [astro-ph.CO]]. 
2] A. Hojjati, L. Pogosian and G. -B. Zhao, JCAP 1004 (2010) 007 |arXiv:0912 .48431 [astro- 
ph.CO]]. 

3] Q. Su, X. He and R. -G. Cai. larXiv: 1204.21461 [astro-ph.CO], 

i] A. Conley, J. Guy, M. Sullivan, N. Regnault, P. Astier, C. Balland, S. Basa and R. G. Carlberg 

et al, Astrophys. J. Suppl. 192 (2011) 1 [arXiv: 1104. 14431 [astro-ph.CO]]. 
5] D. J. Eisenstein et al. [SDSS Collaboration], Astrophys. J. 633 (2005) 560 |astro-ph/0501171| . 
3] A. G. Riess, L. Macri, S. Casertano, H. Lampeitl, H. C. Ferguson, A. V. Filippenko, S. W. Jha 
and W. Li et al, Astrophys. J. 730 (2011) 119 [Erratum-ibid. 732 (2011) 129] |arXiv:1103.2976l 
[astro-ph.CO]]. 

7] D. Stern, R. Jimenez, L. Verde, M. Kamionkowski and S. A. Stanford, JCAP 1002 (2010) 

008 jarXiv:0907.3149l [astro-ph.CO]]. 
i] E. Gaztanaga, A. Cabre and L. Hui, Mon. Not. Roy. Astron. Soc. 399 (2009) 1663 
|arXiv:0807.355T1 [astro-ph]]. 
[29] Y. Wang, Phys. Rev. D 77 (2008) 123525 |arXiv:0803.4295l [astro-ph]]. 

[30] Q. Su, Z. -L. Tuo and R. -G. Cai, Phys. Rev. D 84 (2011) 103519 |arXiv:1109.2846l [astro- 
ph.CO]]. 

[31] A. G. Kim, E. V. Linder, R. Miquel and N. Mostek, Mon. Not. Roy. Astron. Soc. 347 (2004) 

909 [astro-ph/030"4509l . 
[32] A. Ealet et al. [SNAP Collaboration], [astro-ph/0210087] . 

[33] I. Maor, R. Brustein and P. J. Steinhardt, Phys. Rev. Lett. 86 (2001) 6 [Erratum-ibid. 87 

(2001) 049901] |astro-ph/000"7297] . 
[34] D. H. Weinberg, M. J. Mortonson, D. J. Eisenstein, C. Hirata, A. G. Riess and E. Rozo, 

larXiv: 1201 .24341 [astro-ph.CO]. 



11 



